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Abstract 

We consider a two-dimensional electron gas with long range disorder. As- 
suming that time reversal symmetry is broken either by an external magnetic 
field or, as in the case of a delta-correlated random magnetic field, by the 
disorder itself, we derive a supermatrix cr-model. As an intermediate step, we 
provide a microscopic derivation of the ballistic cr-model, and find that certain 
corrections to its usual form may become important. We then integrate out 
degrees of freedom corresponding to short length scales to derive a low-energy 
supermatrix cr-model. We find an extra term in the free energy that couples 
to the correlator of local currents. Use of a proper ultraviolet regularisation 
procedure that preserves gauge invariance indicates that the contribution of 
the extra term seems finally to become irrelevant. Within the scope of our 
analyis, we therefore do not find any deviation of the scaling behaviour of 
the delta-correlated random magnetic field model from that of the conven- 
tional unitary ensemble. We then generalize the discussion to include models 
of even longer-ranged disorder, plus short-range disorder. When the disorder 
is sufficiently long-ranged that the local currents become delta-correlated, a 
new term appears in the free energy that does give rise to logarithmic cor- 
rections to the conductivity. A renormalisation group analysis of the free 
energy yields a scaling form for the diffusion coefficient which contains both 
a positive correction, that represents classical superdiffusion, and a negative 
correction, which is the usual weak localization correction. The fact that both 
corrections are of the same order and opposite sign leads to the interesting 
possibility of a quantum phase transition at weak disorder in two dimensions, 
tuned by the relative strengths of the short and long range disorder. 

PACS numbers: 72.15.Rn,73.20Fz,73.23.-b 
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I. INTRODUCTION 



Since the appearance of the scaling theory of disordered metals PH3|, the theory of 
localization in 2D systems has become well established. Of course, this phenomenon depends 
on symmetries. The addition of spin-orbital impurities leads to the violation of central 
symmetry and hence a situation when the resistivity vanishes at small frequencies. At the 
same time, the delocalizing effect of a magnetic field is apparent in a number of well-known 
phenomena, a natural example being the integer quantum Hall effect in 2D (for a review 
see e.g. Ref. |7j). Another example is provided by the random flux model which describes, 
for example, electron hopping on a bipartite lattice structure with link disorder [@||. In this 
model a tendency towards derealization is displayed as the band center is approached due 
to the existence of a chiral symmetry at the band center P,IH|. 

A question of long-standing interest in this context has been the influence of a random, 
static magnetic field on a two-dimensional electron gas. This is an example of a disordered 
system with a broken time reversal invariance, where the effects of the magnetic field can 
become very complicated. Interest in this problem has been stimulated by its relevance 
to a variety of experimental situations. The solution of the problem might help to clarify 
the behavior of composite fermions for the fractional quantum Hall effect near half-filling 
TT|. In this model, the interacting electron liquid is replaced by a Fermi gas of quasiparti- 



cles, each carrying two flux quanta of a fictitious magnetic field. While the Chern-Simons 
field exactly cancels the external field at half-filling, variations in the electron density due 
to screening of the impurity potential leads to fluctuations in the effective magnetic field 
around the zero value. A similar model applies also in the gauge-field description |12| of 
doped Mott insulators, where slow fluctuations in the gauge field may be approximated by 
a static random field. Further experimental realizations have involved the introduction of a 
random magnet field (with non-zero mean) onto a high-mobility 2D semiconductor through 
an overlayer containing randomly pinned flux vortices |14| or type-I superconducting grains 



I"3| or through a nearby permanent magnet |15 



A theoretical analysis of the random magnetic field (RMF) problem is made technically 
difficult due to the spatially long-ranged nature of the vector potential correlations, which 
result even for short-ranged magnetic field correlations. For example, a straightforward 
application of perturbation theory ]TB|fi7[ ] readily demonstrates the appearance of infrared 
divergencies in the Born approximation for the single-particle Green function. Consequently 
a wide variety of alternative techniques have been applied to this problem. A real-space 
path integral representation has been introduced by Altshuler and Ioffe [[16] and Altshuler 
et al. [0, while the eikonal Jl9| and related paraxial pOf approximations have also been 
employed. If the correlation length, d, of the random magnetic field is sufficiently large that 
d^> 1, 1/kF, where I is the single-particle mean free path and kp is the Fermi momentum, 
then a "classical" regime is reached in which the contribution of classical memory effects 
has been shown to be significant [pl] , p2]| . If in addition the RMF is sufficiently strong that 
d is greater than the average cyclotron radius, then a description of electron transport in 
terms of percolation between a network of 'snake'-states at zero B-field contours becomes 
appropriate |23|-p5|. 

Numerical investigation of the RMF model (equivalent to the random flux model away 
from the band center) has continued without apparent consensus, mainly due to the great 
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difficulty of distinguishing between derealization and localization of states with the very 
large localization length that is typical for two dimensions at large conductances. Conclu- 
sions for this model are divided between localization of all states p6^p8[, the existence of 

and the localization of all states except at precisely zero energy 
It is also of interest that, even for the case of one dimension, recent analytical 
work supports the existence of a metal-insulator transition in the 



a critical region p9|-p4 



H0- 

and numerical 



presence of sufficiently long-ranged disorder correlations. 

In this paper, we study the 2D electron gas with long range disorder. We assume that the 
time reversal symmetry is broken. Our approach applies both to the RMF model and also 
to a model with a long range, potential disorder and a constant magnetic field. The latter 
model is very close to models with short-range disorder with broken time-reversal invariance 
(unitary ensemble) for which the localization of all states is well established. The leading 
order weak localization correction to the conductivity appears, for short range disorder, 
at two- loop order and is negative [j5[ |3l2| ) p[D| . While this result has been also derived from 
conventional diagrammatics | 3~8| , an equivalent and more convenient procedure is provided 
by field theoretical methods based on a nonlinear cr-model. In addition, the use of the o- 
models allows one to prove the existence of the renormalization group (RG) and go beyond 
perturbation theory. 

Aronov et al. [41] have already provided a generalization of the standard field-theoretical 
approach [§,0] to the model of a delta-correlated RMF, in deriving an appropriate form of 
a supersymmetric cr-model. They found that the cr-model was identical to the one derived 
previously for short-range disordered systems with broken time-reversal invariance (the uni- 
tary ensemble) ||40|| , and therefore concluded that all states are localized as for the unitary 
ensemble. In an earlier paper, Zhang and Arovas [42] had proposed the existence of an 
additional term in the cr-model originating from logarithmic interaction of topological den- 
sity or, equivalently, current-carrying edge states around magnetic domains; this term was 
argued to lead to a Kosterlitz-Thouless transition from the localized to extended states with 
power-law correlations. Unfortunately the calculation of Ref. [|2| relies on an incorrect re- 
lation between a correlator of Hall conductivities and the average longitudinal conductivity, 
and the additional term was not found in the analysis of Aronov et al. |[41]| . As concerns 
the models in a constant magnetic field with an arbitrary range of disorder correlations, a 
conventional a-model, with a classical diffusion coefficient appropriate to the correlations, 
has been derived fHJ using similar methods to the work of Aronov et al. 

Even so, the question has so far remained open of whether the new term written by Zhang 
and Arovas |42| can exist in principle and, if so, what would be the possible consequences. 
We examine this possibility more carefully in this paper. First we focus on the case of a 
delta-correlated RMF model. For this model, we provide a derivation of the appropriate 
cr-model 



41 



that is somewhat more complete than the original work of Aronov et al. 
Our derivation demonstrates the appearance of certain technical subtleties that were not 
realised in the original work of Ref. |^T| . This work is an extension of calculations presented 



previously in a short letter Hq] . As a first step we derive a form of the free energy that has 
appeared previously under the name of the ballistic cr-model [pl] , p2| . In fact our derivation 
(also presented in the Letter JI5]) represents the first formally justified derivation of the 
ballistic cr-model, since the long-ranged nature of the disorder provides a small parameter 
(the small ratio of the single-particle and transport lifetimes) to allow a controlled separation 
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of modes and expansion of the free energy, in contrast to previous derivations |5lJ52|. Indeed 



this approach allows us to go further and derive further terms to the usual ballistic cr-model 
that become important for shorter-range disorder. We remark that our derivation of the 
ballistic a-model is valid for long-range potential disorder as well as for vector potential 
disorder, with minor modifications. 

Integration over the degrees of freedom in this model that correspond to length-scales 
shorter than the transport mean free path allows us to derive the final form of the free energy 
that is applicable for the description of low-energy transport. Along with the usual terms of 
the unitary cx-model, we find an extra term which couples to the current-current correlator. 
This term appears also for models of long-range potential disorder, as well as vector potential 
disorder, as long as time-reversal symmetry is completely broken. The new term is similar 
to that written by Zhang and Arovas, except for the appearance of an additional factor II (q) 
that represents the current- current correlator. 

For evaluation of the quantity II(q), it is necessary to employ an ultraviolet regularization 
procedure which should, as shown recently by Gornyi et al. |46j, be formulated carefully in 
order to preserve gauge invariance. The use of this procedure then indicates that the factor 
Il(q) should vanish as q — > 0, as confirmed by diagrammatic arguments. This implies that, 
within the scope of our analysis, the new term remains irrelevant for the delta-correlated 
RMF model in the calculation of corrections to the conductivity that are logarithmic in 
frequency, although it may still lead to higher order corrections known as memory effects 
(we remark that our earlier letter j45| employed a more naive regularisation procedure which 
led to incorrect conclusions). 

As a result we do not find evidence for deviation of the scaling behaviour of the delta- 
correlated random magnetic field model from that of the conventional unitary ensemble. As 
such our conclusions coincide with those of the original work of Aronov et al, although for 
considerably more subtle reasons than were originally realised in that work (which remains 
a partial analysis since they neglect higher harmonics of the Q- matrix). Even so, certain 
questions remain, concerning the role of the massive modes of the field theory, which deserve 
further study before a definitive conclusion can be drawn for this problem. 

We then generalize the argument to consider the possibility of longer-ranged disorder. 
We derive the free energy for a model of short-range disorder plus longer-ranged disorder, 
and showed how the long-range disorder couples to the free energy via a Wess-Zumino term. 
If the long-range disorder is sufficiently long-ranged that correlator of local currents remains 
finite in the limit of q — > 0, then a new term appears in the free energy which is responsible 
for new corrections to the conductivity that are logarithmic in frequency. For this model, 
the new term does takes the same form as the one originally written by Zhang and Arovas 
|f4~2|| (although the latter work does not contain a correct microscopic justification). As 
follows from the preceding discussion, for this term to exist the magnetic field correlations 
must be longer-ranged than delta-correlated: instead, the Fourier transform of the correlator 
(B(r)B(r')) must diverge as 1/q 2 as q — > , where B is the magnetic field. 

At the classical level, this free energy represents particle motion in a field of random veloc- 
ities: similar classical models have already been considered [^8|-p0[], assuming the velocities 
to be delta-correlated. Indeed, Kravtsov et al. |49j and Fisher et al. |5(| have demonstrated, 
via a renormalization group analysis, that particle propagation in this purely classical prob- 
lem (the "diffusion-advection model" ) obeys a "super diffusion" . The electron gas with the 
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long-range disorder represents a quantum generalisation of this model, in which quantum 
coherence phenomena are also present and interplay with the above classical effects. 

By performing a RG analysis of the free energy, we find that two competing contri- 
butions to the scaling form for the diffusion coefficient appear. The positive contribution 
corresponds to classical super diffusion in the presence of long-ranged currents, while the neg- 
ative contribution is the quantum mechanical, weak localization correction for the unitary 
ensemble. Since the classical and quantum contributions to the conductivity are of opposite 
sign and potentially of the same order, we have the interesting possibility of a quantum phase 
transition in two dimensions at weak disorder, tuned by the relative disorder strengths. 

The plan for this paper is as follows. In section [TI], we introduce the model for the delta- 
correlated RMF. In section III, we derive the free-energy functional for this model in terms 
of a supermatrix a-model. This involves deriving the ballistic a-model miscroscopically 
for this problem. We then simplify the free energy by integrating out degrees of freedom 
associated with length-scales shorter than the transport mean free path. The derivation 
is applicable with minor changes to a model with a potential long range disorder and a 
constant magnetic field. In section [IV], we present direct perturbation theory calculations of 
the conductivity and the current- current correlator. In section [VJ, we review classical models 
of diffusion in a field of random velocities, and consider an electron gas in the presence of 
both short-range and longer-ranged disorder. We show how the local current couples to the 
free energy through a Wess-Zumino term, and how averaging over this term leads to a term 



similar to that written by Zhang and Arovas [42|. In section VI we subject the free energy 



to a renormalisation group analysis and derive the scaling form for the diffusion coefficient. 
This demonstrates the appearance of a positive contribution to the scaling form for the 
conductivity in competition with the usual negative, weak localization contribution. Section 
|V11| concludes with a summary and discussion. Some technical details are presented in the 
Appendix. 



II. RANDOM MAGNETIC FIELD MODEL 

As we have mentioned in the Introduction, the considerations below are applicable both 
to a RMF model, and to a model with long range potential impurities and a constant 
magnetic field (to break time- reversal symmetry). While the derivation of the a-model is 
similar for the two types of model, to be specific we carry out calculations for the RMF 
model. We write the Hamiltonian as follows, 




where e and m are the electron charge and mass, c is the velocity of light, ep is the Fermi 
energy and A is the vector potential. We focus below on the case where the magnetic field, 
B = V x A, is delta-correlated, 

(B(r)B(rO> = 2(^)V(r-r'), (2) 

with a strength characterized by the dimensionless parameter 7 (and with (B) = 0), although 
one may as easily consider a more general case with a finite range of the correlations. We 
choose the London gauge for the vector potential A 
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divA =0, A n = 



(3) 



where A n is the component taken at, and perpendicular to, the surface of the sample. The 
sample may be either truly two-dimensional, with the thickness of a single atomic layer, or 
simply have a two-dimensional geometry. While we write explicit formulae for the former 
case, extension to thicker samples is trivial. 

The correlations of the vector potential A corresponding Eqs. (0) and (0) are long-ranged 
and we choose them in the form 

/ >, m 2 r 2 

{A\v)A^v')) = ——V^{v-v'), 

(4) 

The correlator V l i shown here is transversal, that is, 

i 

which corresponds to the transversality of the vector potential A, Eq. (|3|). The parameter 
k 1 is a cutoff that renders finite the otherwise-infinite range of the disorder correlations. 
The 5-correlated fluctuations of the magnetic field correspond to the limit k —>■ 0. The cutoff 
k appears [p^EUD ^ n perturbation theory for the infrared-divergent single-particle lifetime. 
For example, the simple Born approximation yields for the mean free time 

: 2lT dip cos 2 (W2) 
r B ^ = 4e F7 ' 



2vr4sin%/2) + k 2 
2e Fl /K. (5) 



The self-consistent Born approximation (SCBA) displays a weaker divergency, 

^dif 

7T ± ' _/n 27T 

' (6a) 



t scba = z e F7T"scBA / / d£ sin 2 X 



(4snr%/2) + e 2 /e|)(4sin 2 (v9/2) + eje\ + k 2 ) ' 
rsc 1 BA^2e F ( 7 ln(l/K)/7r) 1 / 2 , (6b) 



due to off-shell contributions in the collision integral, Eq. (paj) . The inverse transport time 
r^ 1 remains convergent however in the limit of K — > 0, 

f-271 



tFl l f^-co^cot 2 ^) 
= e^7- (7) 

Although we focus on the above choice of correlator V^(q) in this paper, we emphasize 
again that our method remains valid for an arbitrary form of disorder correlations, with 
minor modifications for scalar rather than vector potential disorder. 
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III. FREE ENERGY FUNCTIONAL 



In this section we derive the free energy functional for the 2D electron gas in the RMF. 
We employ the supersymmetry method which by now has been extensively developed |40 



as a method for the exact evaluation of spectral and wave function properties of metals 
with delta-correlated disorder. In the conventional case the free energy functional takes the 
form of a nonlinear cr-model containing an 8 x 8, position-dependent Q-matrix. In order to 
generalize the method to allow for long-ranged disorder correlations, it becomes necessary 
to employ a Q-matrix that depends on two position variables, rather than one. Following a 
Fourier transformation and a harmonic expansion around the Fermi surface, the Q-matrix 
may be represented as depending on a position variable and an angular harmonic index. 

The degrees of freedom contained in the non-zero harmonics of the Q-matrix are weakly 
massive for a long range disorder and we integrate them out rather than neglect them. Al- 
though the coupling of the non-zero harmonics to the zeroth harmonic is weak, integration 
over the nonzero harmonics requires some care to determine the possible renormalization of 
the effective functional for the zeroth harmonic. One consequence of this renormalization, 
due to the first harmonics only, amounts to the inclusion of two-particle vertex corrections 
to the single-particle lifetime, r, leading to its replacement in the free energy by the trans- 
port relaxation time, r tr . This was the main conclusion of Ref. where arbitrary range 



potential impurities were considered, and of Ref. |41| , devoted to study of the RMF problem. 

We find however a further difficulty associated with the integration over non-zero har- 
monics that was not examined in the previous works (^[^TJ. Namely, additional terms arise 



in the free energy whose evaluation requires a careful integration over higher non-zero har- 
monics, rather than only the first harmonic as considered in the approach of Refs. P^Lffl 



While the contribution of these terms seem finally to be irrelevant to computation of log- 
arithmic corrections to the conductivity, as we discuss below, the proper demonstration of 



this fact is far from trivial and beyond the original analysis of Ref. [41 . 

In order to properly include the potential contribution of the higher harmonics, we derive 
as an intermediate step a free energy that is similar to the "ballistic cr-model" of Muzykantskii 



and Khmelnitskii |51||. The latter model itself represents a generalization of the diffusive 



cr-model to the case of ballistic disorder, where the typical energy scale of Q-fluctuations 
may be as large as the inverse scattering time. In fact we also derive certain additional 
terms to their form of the free energy that become relevant for shorter-range disorder. 

Employing standard methods j|(J, we introduce a supersymmetric ip- field that contains 
eight components, corresponding to fermion/boson, advanced/retarded and time-reversed 
sectors. An averaged product of Green's functions may then be expressed in terms of a 
functional integral weighted by a Lagrangian that is quadratic in the ?/>-field: 

(e./ 2 M)£% 2 (0,r)> 



= -A J ^(r)^(0)^(0)^(r)e- £ P^, 
C = -i J m (Wo + £ A Vr 3 + ^) ^(r)dr, (8) 



where 
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V 2 , (co + i5)A 

and T3 is the Pauli matrix in time-reversal space. 

We now average over the vector potential, A. In doing so we neglect the term in A 2 in 
Eq. (H). This approximation is standard and has been used in previous analytic works on 
the RMF model • Furthermore, for the case of long range potential disorder, for which 

our general method is valid, AMike terms are absent anyway. As for the question of gauge 
invariance, note that our derivation has assumed already a choice of gauge through the form 
of the vector potential correlations (see section [TT]). Nevertheless Appendix |A] contains a 
discussion of how the A 2 term may be handled more carefully and found to be negligible. 

The term linear in A induces after the averaging a term that is quartic in the ^-fields 
and, in contrast to the case of short-range impurities, non-local in position: 



£ int = J (^(r)T 3 V*V<r))V«>-r') ${v') n Vi^{v')) d 2 vd 2 v' . 
Following an integration by parts, L may be rewritten 



C 



int 



^ Q (r)^(r / )r3^^V r , r ,^(r / )^(r)r 3aQ , 



^' = - \ E - ( V r - V r') ( V r " V r') 



(9) 



where k = diag(l, —1) in boson- fermion space. We now decouple C m t via an 8 x 8 matrix 
Q(r, r'). While in the calculation for time reversal invariant systems [50], this step requires 
a careful identification of the slowly varying modes, in the case of the system with the 
broken time reversal invariance, all slow modes are easily identified as corresponding to 
pairs ijj a (r) (r') . Following an integration over the ^-fields, we find the Lagrangian 



C 



--Strln [iH (r)5{r - r') + V r yQ(r,r'] 



+ -Str(Q(r,r')K, r 'Q(r',r; 



drdr', 



(10) 



Here Q = Q„ + iQ±, whereas Q = Q u + Q±_ and Q u (Q±) commutes (anticommutes) with r 3 . 
The Q-matrix satisfies the standard symmetries Q(r,r') = Q(r',r) = KQ^ (r' ,r)K , where 
Q = CQ T C T and C and K defined as in Ref. [gOJ. 

The free energy, Eq. (0), in principle provides an exact description of the system but 
it is not useful in this form as it allows for all possible energy scales for Q-fluctuations. 
As a first step in its simplification we search as usual for the saddle-point value of Q and 
perform a gradient expansion around this value to identify the free energy for low energy 
fluctuations. 
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A. Saddle-point 



To search for the saddle-point to Eq. fllPP, we take Q± = and Q to depend onr-r' only. 
This means that the cooperon "degrees of freedom" are suppressed and that the supermatrix 
Q at the saddle point is homogeneous in space. After Fourier transforming with respect to 
r — r', the saddle-point equation reads Q p = g p , where 

(* - £ P - |A - 2i J ^V^ pi pUQ p ^j g p = -i, (11) 

and £ p = p 2 / (2m) — ep. 

The saddle point of the Lagrangian, Eq. ([TOD, is continuously degenerate and the solution 
of Eq. (|TT|) can be written generally in the form 

Q P = V P A P V P , V P V P = 1 (12) 

where A p is a diagonal matrix depending on p in a complicated way. The analysis simplifies 
however in the limits of k ^> (e^r) -1 (short-range disorder) and K <C (e^r) -1 (long-range 
disorder), the latter case being relevant here, where r is the mean free time to be determined 
from the solution. For either case we may employ the ansatz 

i 

Ap = £ p -e + *A/(2r)' ^ 
where A = diag (1, 1, 1, 1, —1, —1, —1, —1), which leads to the condition 

2/E^- P1 ^ A p^- (14) 



A 
27 



For short-range disorder, k ^> (e^r) -1 , the integral in Eq. ( |i4|) is as usual dominated by 
the on-shell value of p — pi and r takes its simple BA value, as given by Eq. ([5]), which is 
equivalent to the SCBA value in this limit. 

In the opposite limit of long-ranged disorder, k <C (e^r) -1 , which is relevant to the RMF 
model, the off-shell dependence of V p _^ contributes significantly to the integral. In this 
limit, substituting Eq. (|T3| ) into Eq. (|14]) we come to Eq. (^a|), thus obtaining the SCBA 



value of r, given by Eq. (|6b|) . This value of r is less divergent than the simple BA value, 
Eq. (||p. In the following we will denote this SCBA value as simply r. At the same time 
we ensure that we stay within the weak disorder limit, epr 3> 1 (equivalently, 7 <S 1). This 
condition restricts the space of relevant fluctuations in the Q-matrix, as we see in the next 
section. 



B. Fluctuations 

To proceed we expand the free energy in fluctuations of Q about its saddle-point value. 
In the case of short range impurities, fluctuations of the eigenvalues of the supermatrix Q 
are massive with the characteristic energy 1/r and may be neglected. In this case, what 
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remains is to consider massless (in the limit u — > 0) modes j|(J] with energies smaller than 
1/t, which leads to the conventional supermatrix cr-model. 

In the models with long range disorder considered now, however, an additional energy 
scale appears due the difference between the mean free time r and r tr . In this case fluctua- 
tions separate into three types: 

(a) hard massive, with a mass of r -1 , 

(b) soft massive, with a mass of r^ 1 , and 

(c) 'massless', that is, with a mass vanishing in the limit uj — > 0. 

Type (a) fluctuations are associated with the fluctuations of the eigenvalues of Q and the 
Cooperons. The weak disorder condition e^r ^> 1 ensures that these fluctuations are irrele- 
vant and hence we may restrict attention to Q-matrices of the form Q 2 = 1 and [Q, T3] = 0. 
Alternatively, one can use the form of the supermatrices Q given by Eq. ( |T2| ) with [V, T3] =0. 
The effective functional for the remaining (b) and (c) modes will then contain energy scales 
only much smaller than r -1 . 

The above separation of scales may also be expressed in terms of length scales: at 
distances exceeding l tr = VpT^, density relaxation in the classical limit is described by a 
diffusion equation, whereas for distances between I = vft and l tT one should use a Boltzmann 
equation. The possibility to generalize the above classical descriptions to include quantum 
fluctuation phenomena is afforded by the a-model technique. 

To derive a free energy functional describing both soft massive and massless modes we 
proceed as follows. First we introduce the following Fourier transformation for Q(r, r') with 
respect to r - r': 

Q(r, r') = J t^Qp(R) exp(ir 3 p(r - r')), (15) 

where R = (r + r')/2. The appearance of the T3 factor in Eq. ( fTED ensures that Q P (R) 
satisfies the symmetries 

Q P (R) = Q P (R) = KQl(R)K. 

Next we introduce the following parametrization for Q P (R): 

Q p (R) = U(R)qWu(R), 

Qp 0) = K(R)A P K (R) , (16) 

where n = p/[p|. In Eq. ([16]), U (R) is independent of n, and Qp (R) contains only non-zero 
harmonics in n around the Fermi surface. The matrices U and V n obey the symmetries 

uu = 1, U = KWK, 
V n V n = 1, V n = KV^K. { '> 

In what follows we will make expansions in deviations of the supermatrix Qp^ (R) from 
A p . For this purpose the supermatrix Q p °^ (R) can be conveniently parametrized, for exam- 
ple, as 
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with J dnP n = 0, P n = —P n , {-Pn,A} = 0. Both U and V n vary slowly with R, that is, 
on length scales longer than I = vpr. While U(r) represents the massless (c) modes, the 
weakly massive modes (b) are contained in V n . It is the dependence of the supermatrix V n 
on the vector n that leads to the gap in the spectrum of excitations of the (b) modes. 
An important feature of our parametrization for Q p (R) in the form 

Q p (R) = T n (R) A p f n (R) 

is the separation of the rotation matrix T n = UV n into the U and V n factors. This step is 
not only convenient for computation but, more importantly, ensures the preservation of an 
original symmetry of the initial Lagrangian, Eq. ([§]). Specifically, it ensures that the final 
free energy functional is invariant under global rotations U (R) — > UqU (R) with Uq = Uq 1 
independent of coordinates. This invariance follows from the original invariance of the 
Lagrangian, Eq. (|8|), under global rotations ip(r) — ► U ip(r). 

Substitution of the parametrization flTBp into expression (p70|) yields a free energy in 
terms of only the weakly massive (b) modes and the massless (c) modes. Accordingly this 
free energy will describe fluctuations over energy scales only much less than r" 1 . In order 
to investigate the truly low-frequency behavior of transport coefficients, however, we need 
to reduce this free energy even further by integrating over the weakly massive (b) modes. 
This procedure will reduce the free energy to a form in terms of only U(r), representing 
the massless (c) modes. In this way we obtain the final form of the free energy that is 
applicable on energy scales much less than r^ 1 and so appropriate for the description of 
low-energy transport. It is necessary to integrate carefully over the (b) modes, rather than 
simply neglecting them, because, even though their coupling to the (c) modes is weak, it is 
sufficiently complicated to lead on integration to a potentially non-trivial renormalization of 
the bare free energy for the (c) modes. Due to the soft mass of the (b) modes it is sufficient 
to apply a Gaussian approximation in P n ; higher order terms give a small contribution 
provided the inequality ei?T tr 3> 1 (equivalently, 7 1) is fulfilled. 

In principle, if the correlations of the disorder decay sufficiently slowly (for the form of 
the correlations given by Eq. (|]), this applies in the limit of k — > 0), one needs to employ a 
cutoff at large momentum k (short distance). The need for an ultra-violet regularization is 
well-known in the context of the ballistic a-model and originates from the fact that 

the gradient (Liouvillean) operator in the logarithm of the initial Lagrangian, Eq. fllCf ), is 
singular in two dimensions. It turns out that in this model one needs to be careful about 
the precise form of the short- distance cutoff, a point which we discuss further below. 

It is helpful at this point to make contact with the conventional calculation for delta- 



correlated impurities ||40|| . To do so we simply set V n = 1 in Eq. ( |Tq) so that Q becomes a 
function of R only. A straightforward gradient expansion of the Lagrangian fllQf) to second- 
order in the gradients then recovers the conventional a-model. It is useful also to note 
that, prior to the gradient expansion, the Lagrangian (|10D is invariant under a local "gauge" 
transformation 

U(R) -> U(R)h(R), [h(R), A] = (19) 
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(note the term "gauge" is used here in a separate sense from that of the original vector 
potential). As we will see later, this invariance remains conserved also for the case of long 
range disorder and is related to conservation of local currents. 



C. Gradient expansion 



We proceed with the derivation of the free energy functional by expansion of the loga- 
rithm in the Lagrangian ( |T0| ) in low energy terms. Substituting the parametrization, Eq. 
(IHf), into the Lagrangian (|10D , and cycling factors under the supertrace, we find 



C 



Jstrln ( -i£ p + — + A kin + Aoii 
2 \ It 



A 



+i / Str(Q$\R)V^ pi p[p{Q^(R))^ 



(2tt)' 



dR 



d 2 p 

(2tt)» 



(20) 



where 



Am = WF^n(R.)nV R T n (R)r3, 
Aon = 2 J V^ pi p\p{ (f p (R)g pi (R)T p (R) - A P1 ) 
A = -i (oo + id) f n (R)AT n (R), 



d 2 Pi 

(27T) 2 



and T n = UV n . The term Ain describes the kinetic energy, A oll is the collision integral and 
A is the usual frequency term entering the cx-model. Then we perform an expansion of the 
logarithm in the terms Am, Aoii and A- If we keep to first order in all of these terms, we 
find that the contribution of Aoii is (—2) times that of the last term in Eq. (f2~0"|): in this 
way we recover the usual form of the a-model for ballistic disorder: 



J i f i»(n,n')(Q„(r)-Q„,(r)) 



, dndn' 

Ar 3 T n (r)nV r r n (r)— + iu / AQ n (r) — 
Ztt J Ztt 

where v is the density of states and iu(ni, n 2 ) is defined as 



dr, 



w( 



n2 n\n{ 



(21) 



(22) 



so that 



2t B a 



w{ n, n 



A dndn' 

(2tt 



.2 ■ 



1 f dndn' 
2^ = 7 W(n ' n)(1 " nn) (2^- 



(23) 



In the limit k — > the function w (ni, n 2 ) takes a simpler form 
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vf,i/7T7 (1 + nin 2 ) . , , 

ty(ni,n 2 ) = — -, r [k — > 24 

V ' ; 2 (1 -mn 2 ) 



Notice that the free energy F, Eq. (|21|), does not contain any divergencies in the limit 
of k — > 0, since the singularity in w(n, n') as n — > n', Eq. (p4[), is compensated by the 
{Qn — Qn') 2 factor in the collision integral of Eq. (21). This is a reflection of the fact that 



the free energy contains energy scales only less than r since we have eliminated the hard 
massive (a) modes. 

We emphasize that, while the form of the ballistic cx-model in Eq. (ETj) has been presented 
previously |51]]52|| , the derivation here for the case of long-ranged disorder represents the first 
formally justified derivation of this model. This has been made possible by the separation 
of the single-particle and transport lifetimes, r and r tr : this leads to the small parameter 
t/th <C 1 which separates the (a) type fluctuations from the (b) and (c) type fluctuations, 
and allows a gradient expansion of the free energy in a controlled manner. The derivation 
remains valid also for long-range potential disorder, instead of vector potential. In this 
case, Cooperons become operative unless time-reversal symmetry is completely broken by a 
constant magnetic field. 

An interesting point is that, for the RMF model, the ratio T/r tr formally vanishes as 
the cutoff k is taken to zero, at least when the SCBA value of the single-particle lifetime 
r is taken. This suggests superficially that the above derivation of the ballistic a-model 
of Eq. fl2~ip becomes exact for the RMF model in the limit of k — > 0, and in principle 
applicable on length-scales all the way down to the Fermi wavelength. A peculiarity of the 
RMF problem however is that it is not straightfoward to specify the precise value of the 
single-particle lifetime: for example, a gauge-invariant formulation |18| produces a value of 



the single-particle lifetime that is convergent, rather than vanishing, as k — > 0, at least when 
the average magnetic field is nonzero. The precise value of the single-particle lifetime in the 
RMF model that couples to the (a)-mode fluctuations, and hence the value of the shortest 
length-scales down to which the ballistic cr-model is applicable for this model, remains an 
open question. 

A further strength of the derivation of the ballistic cr-model shown here is that, for finite 
values of r/r tr , the expansion of the free energy may be continued to find further terms that 
serve as corrections to the usual form of the ballistic a-model, Eq. (|2T|), but may nevertheless 
become relevant, as we discuss further below. 

Our aim now is to integrate over the P n fields by treating them in a Gaussian approxi- 
mation. To do so we find it convenient to Fourier transform from P n (r) to angular harmonic 
and momentum space via 



J ^ P m ,k exp(i(kr + my)r 3 ) 



d 2 k 

(2vr) 2 



where ip is the polar angle of n. Performing the harmonic expansion on the collision integral, 
we see that, as well as the SCBA scattering lifetime described previously, a whole series of 
lifetimes associated with successive harmonics appears. We define the mth lifetime T^ m > by 

A f d 2 p 



2r( m ) 



y ^I^P- Pl KPi A PiCos(m^), (25) 
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where (p is the angle between p and pi, so that r(°) coincides with the SCBA r and l/r tr 
1/r — 1 / . For example, the collision term in Eq. (ETH becomes 



coll 



7TZ/ 



7tz/ f . . , dndri 



(26) 



While the form of the free energy given by Eq. ( |2"1~D with the collision integral ( |2()| ) is 
sufficient to treat the case of long-ranged disorder correlations (for which r tr ^> r), in the 
limit of short-range disorder certain further terms are also relevant and should be included 
for the integration over the nonzero harmonics. These extra terms represent corrections to 
the usual form of the ballistic cr-model. These terms are found by continuing the expansion 
of the logarithm in Eq. ( [20] ) to include terms of order Al oll , -Akin^coib ^kin^coii an d ^kin- 111 
doing so, we ensure that we include all terms in the free energy that contribute to Gaussian 
order in P n and lead to no more than two gradient operators in each term in the final free 
energy. Collecting the various terms in the gradient expansion together, we come to the free 
energy F = P + P„ + F± + F unit where 

r d 2 k ^ n r n r M (\ i \ 

F ° = ™ / (^ Str [Pn^P-n,-*— {--^) 
H {Pm,\t.P— m— 1,— "I" -fm+l,k-P— m,— k^) j 

/oo 
dr Str [%r 3 A (P m P- m -i + P m+ iP-m) ~ t%A(P m P- m -i ~ P m +xP-m)] , 
m=l 

F± = J dr Str [^r 3 A(P 1 + P_ x ) + i^A^ - P_ x )] , 

Pumt = y J dr Str [D ( VQ) 2 + 2^AQ] , (27) 

and k = k x + ikyT^. The supermatrices $ and Q are defined by 

$ = [7W, Q(r) = U(t)AU(t) (28) 

and ($" L ) are the components of $ that commute (anti-commute) with A. The parameter 
D is equal to D = v 2 f t/2. As an intermediate step we have rescaled P m — > P m r^ m ' jr. 
Eqs. ( |2"T|) now represent the free energy that is applicable for energy scales that are much 
less than r -1 , but may still be as large as r^ 1 . 

In the limit of short-range (delta-correlated) disorder for m ^ 0), we see 

that the P m fields in Eq. ( p7| ) are now infinitely massive. This observation requires the 
presence of the extra terms that we have included in addition to those of usual ballistic 
cr-model, Eq. ([H]). These terms lead to the appearance of the extra (t^/t) factor in P , 
as compared to expression (|26|) . Due to the infinite mass of the P m fields, only the terms 
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in F un it remain in the free energy and the conventional unitary cr-model is recovered in this 
limit, as required. 

As a technical point we remark that our form of the free energy (B7|) differs from the 
ballistic cr-model as written in Eq. in a further sense: the lifetimes t^ 71 ' that appear 
in our expressions are defined through a self-consistent Born approximation, according to 
Eq. (f25[) . This expression allows in general for off-shell contributions from the momenta 
Pi. The collision integral in the model of Eq. flZT|), by contrast, amounts to only a simple 
Born approximation (see Eq. (p3[)) as it includes only on-shell contributions through the 
factor w(n, n'). For the RMF model with a finite n, the distinction is not important for 
sufficiently low angular harmonics, as then the relevant momentum integrals are restricted 
to the energy shell anyway due to the compensation of the singularity in w(n, n') as n — > n' 
by the {Q a — Qn') 2 factor in Eq. (pi]). As an example, for m = 1 this compensation has been 
encountered already in Eq. (|7|) for r tr . For harmonics m of the order of k~ x , however, the 
simple Born approximation becomes unreliable and instead the full momentum dependence 
of Vp_ pi must be accounted for. For such high harmonics, the factor of T^ m ' jr in F becomes 
much larger than unity and hence the contribution of these harmonics is strongly damped. 
The end result of these considerations is that for the RMF model we may apply a cutoff to 
the angular harmonic summation in Eq. ( B7p tom< k~ x . 

D. Integration over nonzero harmonics 

Having derived the free energy functional of Eq. ([27]) for fluctuations at energy scales 
much less than r -1 , the next step is to reduce this form, by an integration over the non-zero 
harmonics P n { r ), to one that is applicable at the lowest energy scales, which are much less 
than r t 7 . In other words, we need to average the terms in the free energy that couple 
the (b) and (c) modes (F tl and F±) with respect to the bare (b) mode free energy (Fq) to 
produce the required renormalization of the bare (c) mode free energy (F un i t ). Equivalently, 
we aim to determine the influence of ballistic electron motion on distances smaller that l tr , 
at which distances the classical limit is described by the Boltzmann equation, on quantum 
interference processes at large distances. 

Since relevant terms in the free energy will contain no more than two gradient operators, 
the relevant contribution from this averaging comes from the second-order cumulant of 
F n + F x : 

F + F„ + F ± ^ -~((F U + F ± ) 2 ) Fo . (29) 

While the contribution from the cross term F n F± vanishes, the terms in F± may be 
eliminated by a simple shift in P±\ which leads to a dressing of the bare diffusion coefficient, 
D , appearing in F unit : D — ► D = VpT tT /2. In terms of perturbation theory, the replacement 
of the bare D by the classical diffusion coefficient D corresponds as usual to the inclusion 
of two-particle vertex corrections. Taking into account only the contributions from F unit 
and Fx corresponds to the calculation of Refs. and gives the conventional unitary 

cr-model with the free energy functional F unit , Eq. fl2~T|), and the classical coefficient D. 

What remains is to evaluate the contribution from the terms in F n . In contrast to the 
contributions from F± and F unit , which involve only the zeroth and first harmonics, the 
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contribution from F n involves correlations between higher harmonics. To calculate this con- 
tribution we need to perform the set of Gaussian integrals that correspond to the integration 
over P m with the free energy functional Fq. This step is no more than an application of 
Wick's theorem in the space of the angular harmonics (see also Aleiner and Larkin |53| for 
a further example) and so requires inversion of the quadratic form in P m contained in the 
bare free energy Fq. 

The full form of the Gaussian integration over the P m fields is made clear if we write out 
the components of P m as 



p m ,= i ~™' k \ ,b. 



-'m,k 



We see from the form of Fq that there are no correlations between negative and positive 
harmonics of I? mi k; defining the column vectors a± = (a±i, a±2, • • ■), b± = (£>±i, 6±2, • • •) and 
similarly for p and a, we come to F = F + + F ~ where 

F o = — [ 7^Rk^±)«±,k + b* ±M L(s ± )b ± M (30) 

Ttr J \*>Tt) 

s± = ± ir 3 k y ). 

Here L is a tridiagonal, semi-infinite matrix with the entries 

r tr r( m ) /l 1 

m,m' T \ T 



L ( S ) ) , = ( Z ~ ZJm) ) ^ m < m ' + is ^rn+l,m' + is*6 m>m / + i, (31) 



where m and ml are positive and 

s = — (/c x 4~ iky j 



After a similar reexpression of F u in terms of the vectors a, b, a and p, one can perform 
the averaging. In the process we use the relations Str$ = Str(A$) = Str(r3$) = that 
follow from the symmetries of U given by Eq. (|TTD: hence only the combination Str(r3A$) 

can enter the final formulae. Integrating over the supermatrices P m in ^-^fV Eq. (|29|), with 

the free energy Fo, Eq. ( p7|) , we obtain a term quadratic in Str(r3A$) with a coefficient 
determined by an integral over k. Indeed we reduce the additional term F c in the free energy 
to the form 

Fc = ~wJ str (r3A$l (r)) str (T3A$J (r ' )} WJ (r ' r ' } rfrrfr '' (32) 



where 



IF (r, r') = / nV 3 'r (r, r'; n, n') T (r', r; n', n) (33) 



dndn' 
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The function T (r, r'; n, n') depends only on the coordinate difference r — r' and may be 
written in the momentum representation as 



r (k; n, n') = V ft' 1 (s)) e ^- m < + (f-x e -«H-m'f) 

V / m,m! V / m,m' 



m,m'>0 



(34) 



where <fi and 0' are the polar angles of the vectors n and n'. An alternative formulation is 
that T (r, r'; n, n') satisfies the following Boltzmann-like equation: 

/ tr nV r r (r, r'; n, n') + [w (n, n") T (r, r'; n", n') — = 5 (r - r') 5 (n - n') , (35) 

J 27T 

where W (n, n') is a function of nn', such that its eigenvalues are equal to the diagonal 

entries of — (t^) X j . As we want to derive a free energy functional 

for supermatrices that vary slowly on the scale of l tv , we assume that the functions $ (r) 
depend more slowly on r than IP- 7 (r, r'). 

Using Eqs. (^) one can check without difficulty that the function U 1 ^ (r,r'), Eq. (j33|). 
satisfies a transversality condition 

V ri IF (r, r') = V^TF (r, r') = 0. (36) 

As the function IP- 7 (r, r') depends only on r - r' we perform the Fourier transformation 
in this variable and using Eq. ([36]) write it in the form 

IF (q) = n (q) - . (37) 

Using Eqs. (|32"D , we write the final form of the a- model as 

F[Q] = yJ St AD(VQ (r)) 2 + 2iuoAQ (r)]dr 

i 2 r 

~ 16 / Str ^ 3A$i ^ Str ^ 3A$i ^ ^ ^' ^ drdr '' ^ 

The invariance of the free energy functional F under the local gauge transformations 
means that it can be written in terms of the supermatrix Q only (without <fr). This may be 
done by employing the transversal form for IP-'(q), as given by Eq. fl37D , in the second term 
of Eq. (|38|). Then we obtain 

f[q] = y J str ( D ( v ^ ( r )) 2 + 2%ujA Q ( r )] rfr ( 39 ) 

11 ^ n ° (9) -M( q) M(- q ) ^ 



256 J q 2 v ^ v ^(2tt) 2 ' 
where 

M (r) = Str (r 3 Q (r) [V X Q (r) , V y Q (r)]) (40) 
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represents the local topological density ||7||. The free energy functional F[Q], Eq. (p9[), has 
a similar form to that written by Zhang and Arovas |^2j (although the latter reference does 
not contain a proper microscopic justification), up to the presence of the factor II (g). 

What is crucial is the behaviour of Ho(q) as q — > 0. If IIo(<7) were to remain finite as 
q — > 0, then the extra term in the free energy, Eq. fl39|), would remain relevant and give 
rise to logarithmic corrections to the conductivity. A vanishing value of Ilo(g) means that 
the new term is not relevant, and gives rise to only higher than logarithmic corrections to 
the conductivity. These statements are made clear in section |V], where we show how the 
quantity IP- 7 (r, r') appears in a direct calculation of the conductivity. Moreover we show in 
section [IV B| that the quantity IP- 7 (r, r') has a direct physical meaning: it is proportional 
the mesoscopic current-current correlation function. Thus the question of the potential 
relevancy of the new term in Eq. (|39| ) amounts to whether or not the Fourier transform of 
the current-current correlator remains nonzero in the limit of q — ► 0. 

In fact the evaluation of the quantity Ilo(g) as q — > is a delicate procedure within the 
framework of the a- model. As shown very recently by Gornyi et al. , a proper evaluation 
requires a careful regularisation procedure that ensures the preservation of gauge invariance. 
Their proposal is to evaluate V (r, r'; n, n') by a semiclassical approximation and to impose 
a short- distance cutoff via a minimum path length. Regularising in this way, the value of 
n (g) turns out to vanish as q — > 0. Their proposal was supplemented by a diagrammatic 
calculation for the current-current correlator which seems to verify the use of regularisation 



procedure; we refer to Ref. [(46[ for further details (we remark that our earlier letter |45 
employed a more naive regularisation procedure which led to incorrect conclusions). 

This demonstrates that, within the scope of our analysis, the new term seems to re- 
main irrelevant in the calculation of corrections to the conductivity that are logarithmic in 
frequency. As a result we do not find evidence for deviation of the localization behaviour 
of the delta-correlated random magnetic field model from that of the conventional unitary 
ensemble. As such our conclusions coincide with those of the original work of Aronov et al., 
although for rather more subtle reasons than were originally realised in that work. 

Even so, certain questions remain which deserve further study before a definitive con- 
clusion can be drawn for this problem. In particular, the role of the massive (a) modes, 
describing fluctuations of the eigenvalues of Q and the Cooperons, has been neglected in 
our analysis and should be considered more carefully. For example, there is some ambiguity 
regarding the precise value of the single-particle lifetime, r, in the RMF model. Although 
the SCBA for the RMF model yields a value of t~ 1 that is divergent as k — > 0, it is likely 
that the true value of r _1 that couples to the massive (a) modes should remain finite; in a 
similar manner, a gauge-invariant formulation (T8| (beyond the SCBA) has been shown to 
produce a finite single-particle inverse lifetime, at least when the average magnetic field is 
nonzero. In the formalism of the field theory, this suggests that an improved saddle-point 
may be found which does not contain any divergencies in the limit of k — > 0. Since the 
single-particle mean free path determines the value of the shortest lengthscales down to 
which the ballistic a-model is applicable for this model, and the evaluation of the new term 
in the free energy depends crucially on the regularization procedure at such short distances, 
further study of the role of the massive (a) modes in this problem seems desirable. 

Although the new term in the free energy, Eq. fl39|), appears not to give rise to logarithmic 
corrections to the conductivity for the delta-correlated RMF, it is still responsible for cor- 
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rections that are higher than logarithmic. Such corrections are known as "memory effects" 



and have been examined recently by Wilke et al. [22]. The derivation of such corrections is 



contained in the calculation of the conductivity in next section. 

We remark that a free energy of the form of Eq. fl39[) arises also for a model of long-range 
potential disorder, as well as vector potential, as long as time-reversal symmetry is broken 
by a constant magnetic field (in which case, a further topological term associated with Hall 
quantization appears). If time-symmetry symmetry is preserved, then the additional term 
in Eq. (|39|) is absent for trivial symmetry reasons. 

The question remains whether the new term in the free energy, Eq. Q5D|), is always 
irrelevant or if such irrelevancy has been specific to our original choice of model. In section 
we show that, indeed, for other models of disorder such a new term can remain relevant, 
and give logarithmic corrections to the conductivity. What is required is disorder of even 
longer range than has been envisaged so far, that leads to a finite value of the current-current 
correlator as q — > 0. This situation leads to the interesting possibility of a quantum phase 
transition at weak disorder. 



IV. PERTURBATION THEORY 

In this section we supplement the derivation of the free energy in the previous sections 



we 



by direct perturbative calculations within the field theory formalism. In section |IV A 
calculate the conductivity. We identify the appearance in this calculation of the quantity 
IT- 7 (r, r'), that appears in the free energy, Eq. (|39|), and show how it leads to higher order 
than logarithmic corrections to the conductivity. We also calculate in section [IV B| the meso- 
scopic current-current correlation function, and find it is directly proportional to IP- 7 (r,r'), 
thus providing a physical interpretation for this quantity. 

While we operate within the framework of a field theory, we emphasize that the calcu- 
lations of this section are equivalent to conventional diagrammatics. We remark also that 



the results for the conductivity of section |IV A| may also be derived via an RG analysis for 



the diffusion coefficient, similar to that of section [VI]: the equivalence of the two approaches 
follows from the Einstein relation between the conductivity and the diffusion coefficient. 



A. Conductivity 

In this section, we demonstrate a direct calculation of the conductivity within the frame- 
work of the field theory. According to standard linear response theory (see e.g. |55j40|] ), the 
response R is defined by the relation of the current density J to the oscillating part of an 
applied vector potential 



lit) 



Following Ref. [40], we write the response as R 1 ^(uj) = J dv'R^ij^ r'), where 



R l Hr, r' 



16n 



E ^ " r ^r> ^( r >;(r)(l " n)u*i V>I(r)$(rO>, 



(41) 



7,5=1 
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and where the averaging is with respect to the Lagrangian C in Eq. @. Here 

K = (-*V r - (e/c) A) jm (42) 

is the velocity operator (A is the static vector potential). On averaging the right-hand side 
of Eq. (pFID , we may pair the ip fields, according to Wick's theorem, in three different ways. 
Ri + R 2 + R3 where 

' 2 Str [k(l - r^ r ,g l \v\ r)] Str [k(l - r 3 )7r^ 22 (r, r') 



This leads to R 

R\ j (r,r 



16tt 



r, r 



i^(r,r' 



167T 



16tt 



Str [*(1 - r 3 )r iig 12 (v', r')k(l - r 3 )7r^ 21 (r, r 
Str [k(l - r 3 )7r^ 12 (r', r)^A:(l + r 3 ) 5 21 (r', r 



(43) 



and where the symbol vr indicates a velocity operator whose derivative acts to the left. The 



Green function is defined by g 



a/3 



r, r 



2(ip a (r)ip^(r')) and, following the introduction of 



the Q-matrix, satisfies (c.f. with Eq. (|TT|)) 



e + e F + 



V 2 

2m 



AU(r,r') 



dr"V r>r „Q(r,r")g(r",r') = -%5(v - r'), 



(44) 



where the function V r y is determined by Eq. flj^). In the limit of short-range disorder, 
t = r tr , the fluctuations of the nonzero harmonics of Q n (r) are strongly suppressed, as 
demonstrated in the derivation of the free energy in section C . In this case, Q(r, r') 
becomes a function of R = (r + r')/2 only and Eq. ( fl"3|) reduces to the relations already found 
in the direct computation of the conductivity for short-range disorder (see Ref . [|0| , Chapter 
8, and Ref. ||56|| ). However, in the limit of long range disorder, when non-zero harmonics are 
not suppressed, an additional term appears that contains the quantity n u (r, r'). 

In the limit of long-ranged disorder, for which r tr ^> r, we may perform an expansion 
of the Green function g(p,R) in the collision term J2 P1 ^p-piPiPi (Qpi ~ Qp)> where the 
momentum p arises from a Fourier transform with respect to r — r': this expansion cor- 
responds directly to the expansion of the free energy in section [III Cj in the quantity ^4 co ii, 
and corresponds to an expansion in powers of r/r tr 1. Therefore to leading order in r/r tr , 
we may take 



<?(p,R) 



e H Qy 

2m 2t 



Since we neglect the hard massive (a) modes (see section p IB[) , we have also that [g, r 3 ] = 
and hence _R 3 = 0. Performing the necessary integrals over the modulus of the momentum, 
we find 



R l l 



K 2 



2e 1 vv%T 



dn 

2^ 



--(Str (A;(l-r 3 )(g n (r)-A) 11 )x 
Str (Ml-T 3 )(Q n (r)-A) 22 ))], 



7T- 



16 



dr 



, dndn' 



i„'3 



-n n 



' (Str [k(l - r 3 ) [Q„(r)(l + v F mV r Q n (r))} 



12 



xA;(l - r 3 ) [Q n ,(r')(l + ^rn'V r ,g n Kr'))] 21 ) 



(45) 
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where now the averaging is with respect to the free energy (|27|). In the limit of t <C r tr , to 
leading order Eq. (^) reduces to include only the following term R, originating from R 2 : 



R 



■i j 



7T- 



[evvF ) 



16 



dr' 



dndn' 



-n n- 



(Str [k(l-r 3 )Q 12 (r)k(l-r 3 )Q^(r' 



(46) 



The response R lJ , Eq. (|46|), is the contribution that was absent in the standard treatment 



of models with short range disorder |56],f40]] . It involves integration over non-zero harmonics 
and may be computed by averaging with respect to the free energy fl2~T|). 



As is consistent with the derivation of the final form of the free energy in section [I I_D 



we employ the parametrization of Eqs. ([16]) and ([18]), and expand each of the Q n matrices 
in Eq. ([16|) up to second order in P n . To reproduce the leading-order quantum correction 
to the conductivity, we also expand in the deviation of the U matrix from unity, although 
we delay this step for convenience. In other words, averaging over P n (r) we derive as the 
first step an expression for the conductivity in a form of an integral over supermatrices Q (r) 
depending on r only. Then, one can compute the integral and obtain the final formulae. 



To zeroth order in P n , the expression (|46|) for R vanishes due to the integrations over n 
and n'. On expanding both of the Q n matrices in Eq. (|46|) to first order in P n , averaging 
over P n , and adding the contribution of zero order in P n from Eq. ( (45]) we find the following 
contributions to R: 



Ri j = 2e 2 uD5 ij 
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(Str(£;(l-T3)(Q-A) n (r))x 



7T 

16 



Str(A;(l-T3)(Q-A) 22 (r))>], 
e 2 u 2 D 2 5 ij J dr'(Str [k(l - r 3 )(QV r .Q) 21 (r')x 

k(l-r 3 )(QV r Q) 12 (v)}) 



(47) 



(48) 



where Q (r) = U (r) AU (r) and the averaging is now with respect to the final form of the 
free energy, Eq. (138] ). We remark that to derive the contribution P& it is necessary to include 
a contraction with the terms in the free energy, Eq. (|27|), that are linear in P n . 

In the derivation of the contributions R a and Rb of Eqs. fl4~7|) and (|48|) , the coupling of 



only the first harmonics, P±i, to the U matrix is involved. As was found in the derivation 



of the free energy in section [D [D| , the role of the first harmonics, beyond that of the higher 
harmonics, is to dress the bare diffusion coefficient, D = with classical vertex 

corrections, so that it is renormalized according to D — > D = v 2 F r tr /2. For this reason we see 
that the expressions for R a and Rb include a renormalized diffusion coefficient, D. Moreover, 
up to this renormalization, they coincide precisely with the corresponding expressions that 
are found for the conductivity in the case of short-range disorder. 

An additional contribution, R c , comes from expanding both of the Q n matrices in Eq. ([46]) 
to second order in P n . On averaging over P n with the free energy functional (|27D , we find 
the new contribution that can be written in the form 



R 1 ' 



32tt 



(el 



tr) 



IF (r,rO (Str [k(l - r 3 )Q 12 (r)k(l - r 3 )Q 21 (r')] 



(49) 
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where the function IP- 7 is given by Eq. (|33|). In contrast to the contributions R a and Rb, the 
term R c does not have a counterpart in the calculation for short-range disorder, since its 
presence relies on the long-range disorder correlations that pertain to the RMF model. In 
contrast to the derivation of the terms R a and Rb, the derivation of R c requires consideration 
of all harmonics, rather than only the zeroth and first harmonics. 

The final step to compute the quantum correction to the conductivity is to expand U 
around unity, according to, for example, 

i - ip^ 1/2 



u -{TTTp) • (50) 

in similarity with Eq. ([78]). In the same time, one should also expand in P the free energy, 
Eq. (p8|). It is important that the new term F c leads to terms starting from P 4 and may 
contribute to the conductivity only in higher orders than considered here. Expanding R a , 
Rb and R c to fourth, sixth and second order in P respectively, and averaging over P, we 
obtain the conductivity 

= a [5^(1 + A a + A b ) + A?] , (51) 

where 

A„ 



A v 

c 



t 2 

—I 2 

128 d ' 


A b = - 


t 2 

I 2 

Q4d d} 


t% 1 


" IF(q) 


d d q 
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(q 2 + A 2 ; 


) (27r) d ' 



and 



d d q 



(q 2 + A 2 ) (2tt) 



(52) 



Here do = 2e 2 vD is the classical conductivity and t = 8/(7ruD). The contributions A a _ c 
originate from the expressions R a - C respectively. 

The contributions A a and A& correspond to the usual weak localization corrections to 
the conductivity for short-range disorder in the unitary ensemble. Since II y (q) — > as 
q — > (see Ref. ||46||), the contribution R c is of higher-order than logarithmic. In fact, 
since n y (q) oc q 2 , R c corresponds to a "memory effect" of the type studied recently in 
Ref. ||22|| , and gives rise to a nonanalytic correction oc \u\ to the conductivity as uj — > 0. The 
contributions A a _ c may all be derived via an RG approach for the diffusion coefficient (see 
section |VTD , as follows by the Einstein relation. 

While we have used the cx-model approach to compute all of the above diagrams, we 
emphasize that in principle it is possible to evaluate all of the above diagrams directly using 
the bare Green functions and without reference to an effective Lagrangian. Such a calculation 
would be useful as a check for the regularisation procedure used for the ballistic cx-model. 
The calculation is not simple as, for example, it becomes necessary to take into account 
a proper dressing of the Hikami boxes and vertices with impurity lines H. In addition 
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extra, mutually cancelling diagrams would appear that correspond to the hard massive 
modes of the field theory. The diagrammatic arguments of Gornyi et al. [[|6|] do not fulfill 
the purpose of such a check: although their calculation of the current-current correlation 
function seems valid (see the next subsection), their calculation of the conductivity relies 
on a relation between the conductivity and the current-current correlation function which 
contains the quantity r. As we mentioned above, this quantity formally vanishes within the 
SCBA for the delta-correlated RMF model. Thus their arguments are questionable for the 
delta-correlated RMF model, although they may be applicable when the RMF is a weak 
perturbation to sufficiently strong, short-range disorder, in which case a well-defined r does 
exist. Instead, what is required is a direct diagrammatic calculation of the first logarithmic 
correction to the conductivity, a task which is not simple even for short-range disorder and 
has not been attempted so far for the RMF model. 



B. Mesoscopic Current Correlations 

In this section we continue the perturbative analysis to a calculation of the current- 
current correlation function for long-ranged disorder, such as a delta-correlated RMF. This 
enables us to provide a physical interpretation of the quantity IP- 7 (r,r'), as being directly 
proportional to the current- current correlation function. 

We employ the definition 

J\v) = \im(K - K,)(Gf(v,r') - G?(r,v% (53) 



for the local current, where %\ is the velocity operator, Eq. Q4"2]). The current-current 
correlator 



I^(r-r') = (J'(r)J^r')) imp (54) 



is then given by 



e 2 



^'(r -r'H-;— lim ,((l-r 3 ) 77 «-7r^(r)^(s)x (55) 

OZTT s^r,s'^r' 

y<(l-r 3 )ss(K'-<^Kr')^'))c 



where a summation over 7 and 5 is implied. The angle brackets in Eq. ( 541) stand for 
averaging over impurities, while the angle brackets in Eq. fl5"5] ) stand for averaging with the 
Lagrangian C, Eq. (H). 

There are now three possible ways of making pairings of the ip and ip fields. Note that 
the calculation is following very similar lines to those of the calculation of the conductivity 



in section [IV A| . As in that case, only one of the pairings need be considered at leading 
order, where the first ip and ip are paired, and the second ip and ip are paired. Following 
Fourier transformations with respect to r — s and r' — s', we find 

J^(r-r') = --^ lim Str(*(l - r 3 )(< - ^"(r, s)) x (56) 

xStr(Ml-r 3 )(^-^)/ 2 (r',s')), 
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where g(r, s) is defined as in section |I V A| , Eq. (0). Following Fourier transformations with 



respect to r — s and r' — s', integrating over the modulus of the momentum, and keeping 
the leading-order term in the limit of r <C r tr , we arrive at 

J*(r - r') = ~^v F f j ^n^Str^l - r 3 )Q n (r)(l + A)) x 

xStr(fc(l-r 3 )g n ,(r')(l-A)). (57) 

Eq. ( |57D represents the required formula from which the current-current correlation function 
may be computed (compare with Eq. ([46]) for the conductivity). As for the conductivity, 
we employ the parametrization of Eqs. flIB| ) and (|T8"D, and expand each of the Q n matrices 
in Eq. (|57|) up to second-order in P n . If we expand both of the Q n terms in Eq. (E7h to 



second-order in P n , we find the contribution 

,2 r AnJni 



J?(r - r') = ^ F ) 2 | *^ n W'Str(*(l - r 3 )P n (r) 2 (l + A)) x 

xStr(A;(l-r3)P n ,(r / ) 2 (l-A)). (58) 
Performing the contractions in Eq. (|5"8Q , we find 



/o ? '(r-r') = |^/ t 2 r n^(r ) rO (59) 

where the function IP- 7 (r, r') is given by Eq. fl33l) . We see from Eq. ( [59"D that the correlator 
of mesoscopic currents 1$ depends on the same function IP- 7 that enters the action, Eq. ( P2"| ) . 
and the conductivity, Eq. (|49|) . 

Strictly speaking, the function II 4 - 7 that appears in the free energy represents not precisely 
the current-current correlation function, but the contribution of only the nonzero harmonics 
to this quantity: this is reflected by the fact that P- 7 contains further contributions, be- 
yond Iq , which are zeroth order in P n , but second-order in the generator P of U (see the 
parametrization of Eq. (|50|) ). In a similar way, the function V (k; n, n') of Section [III D| is 



defined in Eq. (^) through a projection onto only nonzero harmonics m > 0. An alternative 
definition of T (k; n, n') may be written which treats the zeroth and nonzero harmonics on 
an equal footing, as in Ref. ||46|| , although the difference becomes inessential for momenta 



k ^> ; the deviation of the function IP- 7 from the true value of the current- current cor- 
relation function is therefore minor. The separation of zeroth and nonzero harmonics is in 
fact inessential in the computation of correlation functions, as performed in this section and 



in Ref. [46[ ; however, it is unavoidable for the derivation of a free energy in terms of a local 
Q(r) matrix that represents only the zeroth harmonics, as has been performed in section 

ra. 



V. CLASSICAL SUPERDIFFUSION 

We now generalise the discussion to include the possibility of longer-ranged disorder 
than has been represented so far by the delta-correlated RMF. The question that arises 
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from the preceding analysis is whether a new term, of the type represented in the free 
energy of Eq. (|39"1), may exist and still remain relevant to give new, logarithmic corrections 
to the conductivity. As has been demonstrated in section fj fD| the potential relevancy of 
such a term is related to the question of whether the current- current correlator remains 
nonzero in the limit of q — > 0. Although such a situation does not arise for the model of a 
delta-correlated RMF, it is not ruled out for models of longer ranged disorder. 

In fact, the problem of diffusion in a field of currents that are correlated in precisely 
such a way, ie. so that the current-current correlator remains nonzero as q — > 0, has already 
been examined on a classical level Pffi-|oTi|]. These authors considered a model of classical 
2D diffusion in a random stationary velocity field, as governed by a classical Fokker-Planck 
equation 

[d/dt + V (v-A)V)] P (r,t) = 0. (60) 

Here P (r, t) is the distribution function of a randomly walking particle, where D is the 
classical diffusion coefficient. The randomness of the velocity v makes the problem in 2D 
quite non-trivial. One can imagine different forms of the velocity-velocity correlations and 
so, different models leading to different types of diffusive behavior were considered. In one 
of the models (Model II in classification of Ref. the following correlation was assumed: 

(v i (r)vi(r'))= l0 K»(r-r'), (61) 

where the function (r — r') is defined as 

K ij (r - r) = 5 lj 5 (r - r) - — V* Vj, In |r - r'|, (62) 

2n 

so that in momentum space, it takes the form 

K ij (q) = & - (63) 

The model of Eqs. ( |5"0"D and (ET]) may describe, for example, Brownian motion of a particle 
in an incompressible liquid with random stationary flow. Although roughly speaking the 
function describes short-range correlations of the liquid, more precisely long ranged 
correlations are present in Eq. fl6~T|) due to the incompressibility. 

In Refs. [f|^,[5(J a Green function of Eq. ( pOf) was introduced and expressed as a functional 
integral over a vector field (in Ref. [jiJj a replica formulation was used). In this way 
calculation of the diffusion coefficient D is reduced to a study of <p 4 theory. In 2D the 
corrections to the bare diffusion coefficient are logarithmic, and hence a renormalization 



group treatment may be employed, up to the two loop approximation in Ref. |j49| . The 
dependence of the diffusion coefficient D on frequency to was obtained in the form 



D (u) = D n 



1 + 9o In 



D k 2 



UJ 



1/2 



(64) 



where k is an ultraviolet cutoff and go = Jo/ (47tDq). The dependence D (u) given by Eq. 
(0) is classified as superdiffusion because the diffusion coefficient grows as the frequency u 
vanishes. 
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With this background in mind, we return to the problem of electron motion in the pres- 
ence of long-ranged disorder. To make closer contact with the diffusion-advection problem, 
we in fact consider electron motion in the presence of both short range and long range dis- 
order. We assume that at short distances the electron motion is determined by the short 
range potential and is diffusive. The long range potential, which may be scalar or vector 
potential or a combination of both, influences the diffusion at large distances and may be 
considered random force. 

Remarkably, we find a direct relation of this electron model with the diffusion-advection 
model, described by Eq. (|5TiD, as far as certain terms in the solution for the diffusion coeffi- 
cient are concerned (see Eq. (|64])). In the electron model, we find two different contributions, 
of opposite sign, in the solution for the diffusion coefficient. The positive part of the con- 
tribution corresponds precisely to the term appearing in Eq. ( ]6^ ) and hence corresponds 
to classical superdiffusion. It derives from a new term in the free energy, which is of the 
precisely the form written in Eq. (^), although with a finite value of the current-current 
correlator (U (q)) as q — > 0; the term therefore is of the same form as that written originally 



by Zhang and Arovas ||42j| . The electron model differs however from the diffusion-advection 
model in that it is quantum mechanical: hence an additional negative contribution appears 
in the solution for the diffusion coefficient. This negative contribution is the usual unitary 
weak localization correction due to quantum interference, and competes with the above clas- 
sical correction. This competition leads to the interesting possibility of a quantum phase 
transition at weak disorder, tuned by the relative disorder strengths. 

In order to derive a field theory for the electron model we apply the formalism developed 
in the preceding sections. We proceed by averaging first over the short-range impurities, 
but not over the long range potential. A local, quartic term appears in the Lagrangian in 
the usual way, which may be decoupled with a local Q(r) matrix. We find the Lagrangian 

C = / B Strh (L ( P " Z Ml) ~ ' T ^T +V(r)-e F + f + !«(,)) 

dr, (65) 



+ ^StrQ 2 (r) 
8r 



where r is the mean free time for scattering from the short-range impurities (corresponding 
to Dq in Eq. fl60|)). The vector potential A (r) and the long range potential V (r), both of 
which may be random, are not specified now. 

Next we write Q (r) = U (r) AU (r), cycle the U and U factors under the supertrace and 
perform a gradient expansion in U"VUt^. To first order in the expansion, we find the term 



s/ dr H str ( ew W T *w(p-j A ))' < 66 > 

where the matrix function g p is defined as 



6F= 

2m 



— (p- -A(r) -ir 3 V, 
im \ c 



2 i A 



£ P (r) 



Note that in the absence of the long range potentials, the term SF would vanish by symmetry 
on integration over p. We notice further that the term 5F may be expressed in terms of 
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the local classical current density J, which has been defined in Eq. 
transformations, J may be reexpressed in terms of g as 



). Following standard 



J l (r) 



d 2 p 

2vrm J (2tt) 2 



where there is no summation over a, which allows us to rewrite Eq. (|66|) in the form 



5F 



71 



— / dr J*(r)Str(Ar 3 $ l (r)). 
2e 



(67) 



where $(r) = U (r) VZ7(r). Eq. (p7|) clarifies immediately the physical meaning of the 
supermatrix showing that it plays the role of an effective vector potential. If the par- 
ticle number is conserved, the current is transversal (divJ = 0), which means that only the 
transversal part of $ may enter the final free energy functional. 

Derivation of the remaining terms in the free energy continues as usual: the gradient 
expansion is continued up to terms second-order in <E», and first order in uj. If we include, for 
the sake of generality, the possibility of a nonzero average component (Bq) to the magnetic 
field, we come to the free energy 



F[U] = Y J ^ r Str(D[$, A] 2 + 2D xy T 3 A$^ + 2iuKQ) 
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+ - / rfrJ i (r)Str(Ar 3 <K i (r)), 



(68) 



where D = and D xy = cBqtD jmc (and we have assumed that the magnetic field Bq 

is classically weak, ie. eBor/mc < 1). The free energy may then by rewritten in terms of 
the Q-matrix, the last term in Eq. (|68|) by means of a reexpression as a Wess-Zumino term 
IHjg (using the identity divJ = 0): 



F[Q] = y J dr Str(D(Vg) 2 - D xy r 3 Q[V x Q, V y Q] + 2iu\Q) 
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+ — I dr J J duStr ( t 3 Q 



du 



(69) 



The free energy ( p9f) represents the most general free energy for electron diffusion in the 
presence of an external magnetic field and nonzero local currents. 

The relation of the free energy ( |69D to the classical models of diffusion in a field of random 



velocities p8|-|50|] is made clear if we consider the correponding classical diffusion propagator. 
To derive this propagator, we consider Q to be close to A and expand to quadratic order 
around A, via Q = A(l + 2iP — 2P 2 + ...). The free energy reduces to 



dr Str P 



-DV 2 -iu Ar 3 JV 

ev 



P 



Comparison of the square brackets in Eq. ( |70|) with the Fokker-Planck equation (]6 
use of the classical relation between the current and velocity, 



(70) 



and 
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J = e^v, 



shows that the free energy (|69|) coincides with the diffusion-advection model on the classical 
level. The free energy ( |B5D however contains additional physics to the diffusion-advection 
model, since it describes also quantum interference phenomena as represented by interaction 
of the diffusion modes and contained in higher powers of P. 

We now average over the random potentials. As we want in this section only to make the 
connection with Refs. p8|-^0|] , we do not require a microscopic specification of the long-range 



disorder correlations, but may simply assume that the current densities J (r) are random 
and delta-correlated, so that 

(J 4 (r)^(r / )) = e 2 G ^'(r-r / ), (71) 



where K l ° (r — r') is defined as in Eq. (|62|). In particular, lf y (q), as given by Eq. (|63"D, 
remains nonzero in the limit of q — > and takes the required transversal form. The term in 
the free energy, Eq. (BBf), that couples to the currents becomes 



7T 2 f 

5F -> -— / drdr'{J'(r) J J (r'))Str(r 3 A$ l (r))Str(r 3 A<^ (r')) (72) 
8e z J 

In this way, we arrive at the free energy 

F[Q] = Y J Str [^(VQ(r)) 2 - D xy r 3 Q[V x Q, V,Q](r) + 2icoAQ(r)]dr 

~tL I drdr ' Ki ° ( r - r ') str to) str ( T3A ^ ( r ')) ' ( 73 ) 

where the dimensionless coefficient (3 is defined by 

13 = 2tt 3 G . (74) 
An alternative reexpression of the free energy of Eq. (|73[) is 



F[Q] = yJ Stt[D(VQ (r)) 2 - D xy r 3 Q[V x Q, V„Q](r) + 2zo;AQ (r)]dr 

-J^Y J M{r)\n\r- r'\M (r') drdr', (75) 

where M(r) is defined as in Eq. fl4Tj|). 

So far we have specified the long-range disorder only through the resultant current- 
current correlations, requiring that the correlator (|7Tf ) remains nonzero in the limit of q — >• 0. 
Although such a formulation has the advantage of generality, it is useful to identify its 
relation to specific miscroscopic models. It is clear from the preceding sections that a delta- 
correlated RMF does not provide the required form of current-current correlations. Instead, 
the magnetic field correlations need to be of longer range: the Fourier transform of the 
correlator (B(r)B(r')) must diverge as 1/g 2 as q — > 0. 

We see that the final term of Eq. (|75|) is of the same form as that written by Zhang and 



Arovas |4"2"||, although the latter work does not contain a correct miscroscopic derivation; 
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specifically, a proper treatment along the lines of Ref . , since they deal with short-range 
magnetic field correlations, would result in a replacement of the function K^(q) appearing 
in the free energy of Eq. (|73| ) by q 2 K l ^(q). Indeed we may apply the methods of this 
section to a model of short-range disorder plus a delta-correlated RMF, to find just such a 
replacement in the free energy of Eq. ([T3|). It is straightforward then to verify that this free 
energy reproduces precisely the results of Wilke et al. [22] for the memory-effect correction 
to the conductivity in this model. 

In contrast with the final term of the free energy, Eq. (|39|), for the delta-correlated RMF 
problem, here the new term does contribute logarithmic corrections to the conductivity. In 
the next section we sum all such logarithmic corrections by subjecting the free energy of 
Eq. ( |75|) to an RG analysis. 



VI. RENORMALIZATION GROUP ANALYSIS 

In this section we subject the free energy of Eqs. (|73|) and (|75|) to a renormalisation group 
analysis. In this way we find that the new term in the free energy (the final term of Eq. (|75| ) ) 
leads to a new contribution in the scaling form of the diffusion coefficient with respect to 
frequency. The new contribution is positive in sign, and hence competes with the usual 
negative contribution that represents weak localization. Moreover the new contribution has 
precisely the same form as that which appears in scaling form in the problem of classical 
superdiffusion p^-jSP], according to Eq. (|64]). 



Our approach is entirely analogous to the standard procedure applied in the case of 
unitary disorder p9|,|3|,^0| . The only differences are caused by the presence of the additional 



term in the free energy (|73|). The latter term leads to an extra 'effective charge', /3, whose 



flow is coupled to that of the usual charge t = 8/(iri , D). In the standard case, the first order 
(one-loop) contribution in the conventional unitary a-model vanishes and the localizing 
behavior originates from the two- loop diagrams that are the next order in t -1 . For the free 
energy (|73|), the /3-term contributes to the effective charge t already in the first order and 
therefore this contribution can be larger than the conventional one. While the calculation of 
the conventional two-loop diagrams is rather involved and needs a special regularization (a 
dimensional regularization is usually used), the contribution of the /3-term is much simpler 
and does not need any such regularization. 

Following a standard procedure (see e.g. Ref. flOfl ) we separate the unitary supermatrix 
U(r) into slow U(r) and fast Uo(r) parts: 

U(r) = U{t)U {t). (76) 

The fast part Uq (r) contains in the momentum representation momenta in the interval 
A < k < k , where k is the upper cutoff momentum. To guarantee that the rotational 
symmetry is not violated after this integration we impose an infrared cutoff in the momentum 
integrals by adding the following term in the free energy 

F re 9 = -— J Str (AQ) dr, 

where A 2 ^> uj = ujjD. 
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Substituting the parametrization, Eq. (|76|) , into the free energy (|73|), we find F = 

F l + FP, 

F l = \J rfrStr [(Vg ) 2 + 2[Q , VQ ]$ + [Qo, $] 2 (77) 

+2iuUAUQ , 

F/3 = ~T&^ / rfrrfr '^ ( r " r ') [ Str ( r 3^o (r) (r))Str(r 3 Q (r') % (r')) 

+Str(r 3 A$ 0i (r))Str(T 3 A$oi (O) + 2 Str(r 3 Q (r) (r))Str(r 3 A$ 0j (r')) 

where <fr = UVU, <E»o = UqVUq and Qo = UqAUq. Notice we do not include any contribution 
from the topological term (that couples to D xy ) in the free energy of Eq. (|73|), since as usual 
this term does not contribute to the RG flow equations to any perturbative order. 

Our task is now to integrate exp (-F) over U using perturbation theory and obtain a 
new free energy functional F containing second powers of $ Q (r) . We parametrize Uq by 

expand the free energy to the required order in P and average over P to obtain the required 
renormalized free energy. 

It is important to mention that the bare free energy, which is quadratic in P and with 
which the subsequent averaging will be done, originates from F*. As concerns the contribu- 
tion from F 13 , its expansion in P begins with P 4 terms. Therefore the new term does not 
modify the diffusion equation itself but changes interaction between the diffusion modes. 

Terms coming to the renormalized functional from F* are well known and we simply 
use the results written in Ref. ||40|| . It is important that one does not obtain from F* any 
contribution to the renormalized /3-term. As concerns we write F 13 = P-f + P' 3 , where 

if = JL [ drdr'K ij (r - r') {Str(r 3 AP (r) $ 4 (r))Str(r 3 AP (r') ^ (r')) (79) 
Air J 

+Str(r 3 A$, (r))Str(r 3 AP 2 (r') ^ (r')) + Str(r 3 A$ l (r))Str(r 3 AV^P (r') P (r'))}, 
F 13 = -j?r- J dvdv'K ij {v -r , )Str(r 3 A$ J (r))Str(r 3 A<l j (r')). 



Then we integrate over P using the contraction rule j|0 



< Str (PMx) Str (PM 2 ) > = g ^2 + A2) Str (M X M 2 ) , (80) 

and < PP >= 0, where the supermatrices Mi and M 2 anticommute with A and are self- 
conjugate. 

We see that in the first order approximation, there is no contribution to the renormalized 
/3-term coming from F&. So, we conclude that in the one-loop order the /3-term is not 
renormalized and the coefficient (3 keeps its bare value. 
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The first order contribution to the diffusion (first) term in the action, Eq. ( j38|) , comes 
from only the first term in F± in Eq. (|79|) , while the contributions from the second and 
the third terms in Eq. ( |79|) vanish. The contribution from the first term in F± contains the 
integral 



(ft 
32?r 



Str ($f (r) ®f (r)) dr 



d 2 k 

(2^)2 



8 ij 



k 2 J k 2 + A 2 



(81) 



After performing the above integral, the renormalized free energy F can be written in 
d = 2 + e dimensions (for small e) as 



str(vg) s 
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2icDStr(AQ) 



16tt 



drdr'K ij (r - r / )Str(r3A$ ^ (r))Str(r3A<l j (r , )). (82) 



The second term in the figure brackets in Eq. (|82|) is the conventional two-loop contribu- 
tion of the unitary a-model. To write it properly one may use a dimensional regularization, 
which is the reason this term is written for e ^ 0. At the same time, there are no ambiguities 
in the calculation of the third term in the figure brackets and it is written for e 1 already 
in 2D. 

The Gell-Man-Low function for 



t = 

becomes, after continuation to 2 dimensions, 

dt 1-3 
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(3 = const 



t 6 ((3-l) + 0(t A ). 



(83) 



The solution to the flow equation, Eq. (|83D may be found by integrating over A up to 
the ultraviolet cutoff of 1 /Z tr , and we obtain 



t(u) = t c 



l + -G0-l)i|jln 



UT tl 



-1/2 



(84) 



where to is the bare value of t proportional to the classical resistivity. Rewriting Eq. (R|j) in 
terms of the diffusion coefficient and reintroducing G (see Eq. (|7lD ), we find 
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In 
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(85) 



The classical relation J = ev\ between the current density J and the velocity of a single 
particle, v, leads to G = v 2 ^$. The coefficient of the positive correction to the diffusion 
coefficient can then be reexpressed as 
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G 



where go — 7o/4vrZ^Q. Comparison with the scaling equation (p^ ) shows that this coefficient 
coincides precisely with that found previously [f|^,[5(| i n the scaling form for the problem of 
classical superdiffusion. 

This correspondence provides immediately an classical physical interpretation for the 
positive part of the correction to the diffusion coefficient in Eq. fl85|). For the flow to be 
correlated according to Eqs. fl6T|) and (|7T|), random stationary flow 'cycles' should exist in 
the system. A tendency towards ballistic motion along the cycles of local current (see Fig.|l]) 
leads to the positive, or super diffusive, correction to the diffusion coefficient. 
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The difference of the electron model we consider here with the models of classical su- 
perdiffusion [^8-50] is that for the electron model quantum interference processes are also 



important and give rise to negative corrections to the diffusion coefficient. This is the nega- 
tive contribution in the scaling form of Eq.(R5|), and it coincides with the usual contribution 



known for the unitary ensemble for short-range disorder |39], 

The presence of both classical and quantum corrections to the conductivity, which are 
both logarithmic and therefore in competition with each other, provides the interesting 
possibility of a quantum phase transition in two dimensions at weak disorder (at (3=1), 
tuned by the relative strengths of the short-range and long-range disorder. 

In the case when the positive contribution in the scaling form for the diffusion coefficient 
dominates over the negative contribution {(3 > 1), Eq. (|4]) can be used for an arbitrarily 
low frequencies as soon as to <S 1- This is because the effective charge becomes even smaller 
when u — > and one may safely neglect higher order terms in Eq. fl83|) . This type of behavior 
is similar to that known in the symplectic symmetry class that describes short-range, spin- 
orbit impurities M, although the final forms of the scaling functions in the two systems 
are different. In such cases, the type of the renormalization group equations is known as 
a 'zero-charge' situation. The possibility to use Eq. ( |3lD down to zero frequencies means 
that this represents the complete solution. This contrasts with, for example, the localization 
behavior for the orthogonal ensemble where the effective charge (resistivity) t grows when 
decreasing the frequency and the solution of the RG equations loses its applicability as i 
becomes of the order of unity. 

The result for the scaling of the diffusion coefficient with frequency, Eqs. (|85|), may be 
obtained also perturbatively order-by-order, by calculating the conductivity directly, as in 



section [IV A| . In such a calculation a correction, da 1 ^, appears that is similar to R l J of 



Eq.([49]): it is proportional to the current- current correlator and has no counterpart in the 
conventional calculation |56"|J40|| for short-range disorder: 



So* = ^ Jd(r - r')Str(fc(l - r 3 )Q 12 (r')k(l - r 3 )Q 21 (r))< J l '(r)^(r')). (86) 
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Expanding each Q in the right-hand side of Eq. fl86|) to first-order in P, where Q = A(l + 
iP)(l—iP)~ 1 , and averaging over P leads to a correction to the conductivity that corresponds 
to the positive contribution in the scaling form for the diffusion coefficient of Eq. (|85|), after 
a perturbative expansion. The correction to the conductivity corresponding to the negative 
contribution in Eq. fl35|) can also be derived diagrammatically in the usual way |56| , f4*0 



The equivalence of the direct calculation of the conductivity with the RG approach for the 
diffusion coefficient follows naturally from the Einstein relation, which must hold for the 
electron problem. There does not seem to be a connection between the mobility studied 
in Refs. |IH|-^D] and the conductivity of the electron gas with long-range disorder, and its 
dependence on the frequency differs from that for the diffusion coefficient. 

In summary, in this section we have subjected the free energy (|75|) to a RG procedure 
to derive to scaling equation for the diffusion coefficient. As showed in section M, the free 
energy ( f75|) represents a quantum generalisation of the diffusion-advection problem and 
applies when the current-current correlator due to long-range disorder remains finite in the 
limit of q — > 0. The scaling form for the diffusion coefficient contains a positive contribution, 
which corresponds to classical superdiffusion, and a negative contribution, which is the usual 
weak localization correction for the unitary ensemble. Since they are of opposite sign, the two 
contributions may compete and this leads to the possibility of a quantum phase transition 
in two dimensions at weak disorder, tuned by the relative disorder strengths. 



VII. SUMMARY AND DISCUSSION 

In this paper, we have studied the 2D electron gas with long range disorder and broken 
time-reversal symmetry. We have provided a careful derivation of the free energy for this 
system. We have shown how the the ballistic cx-model may be derived in a controlled manner 
for this problem, due to the presence of the small parameter r/r tr , thus providing the first 
formal derivation of this model. We also derived certain additional terms which are beyond 
the usual form of the ballistic cr-model but become important for shorter-range disorder. 
Integration over nonzero harmonics of the Q-matrix leads to a free energy with with an 
extra term beyond those of the conventional unitary a-model. The evaluation of the extra 
term requires the use of a proper regularisation procedure of the ballistic cr-model that 
ensures the preservation of gauge invariance and that has been verified by diagrammatic 
means. 

This demonstrates that, within the scope of our analysis of the delta-correlated RMF 
model, the extra term in the free energy seems to remain irrelevant in the calculation of 
corrections to the conductivity that are logarithmic in frequency. As a result we do not find 
evidence for deviation of the localization behaviour of the delta-correlated random magnetic 
field model from that of the conventional unitary ensemble. As such our conclusions coincide 
with those of the original work of Aronov et at, although for considerably more subtle 
reasons than were originally realised in that work (which remains a partial analysis due to 
their neglect of the higher harmonics of the Q-matrix). 

Certain questions for this model remain regarding the role of the massive (a) modes, 
which describe fluctuations of the eigenvalues of Q and the Cooperons, and which have been 
neglected in our analysis. For example, the precise value of the single-particle lifetime that 
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couples to the massive modes (a) remains to be determined, and is likely to remain finite as 
the infrared cutoff k is taken to zero. This suggests that an improved saddle-point may be 
found which does not contain any divergencies in the limit of k — > 0. A proper treatment of 
the massive (a) modes would allow the application of the ballistic cr-model down to distances 
of the order of the single-particle mean free path, which is desirable since the evaluation of 
the new term in the free energy depends crucially on the the regularization procedure at 
such short length-scales. 

The derivation of the free energy was supplemented by direct perturbative calculations 
within the field theory formalism. We calculated the conductivity directly and showed the 
appearance of an additional correction, that is analogous to the new term in the free energy 
and is higher order than logarithmic. We also calculated the current-current correlation 
function, thus providing a physical interpretation for the coefficient of the new term in the 
free energy. 

We then generalised the discussion to include the possibility of longer ranged disorder. 
We derived the free energy for a model containing both short-range disorder and longer- 
ranged disorder, and showed how the long-range disorder couples to the free energy via a 
Wess-Zumino term. If the long-range disorder is sufficiently long-ranged that the correlator 
of local currents remains finite in the limit of q — > 0, then a new term appears in the free 
energy which is responsible for new corrections to the conductivity that are logarithmic in 
frequency. We have demonstrated how this free energy may be described as a quantum 
generalisation of the diffusion-advection model. By performing a RG analysis of the free 
energy, we have shown how two competing contributions to the scaling form for the diffusion 
coefficient appear. The positive contribution corresponds to classical super diffusion in the 
presence of long-ranged currents, while the negative contribution is the conventional, weak 
localization correction for the unitary ensemble. The competition between the classical and 
quantum contributions to the conductivity leads to the interesting possibility of a quantum 
phase transition in two dimensions at weak disorder, tuned by the relative disorder strengths. 

We remark that a direct diagrammatic calculation of the leading logarithmic correction to 
the conductivity would be a useful check of the predictions of the field theory for the delta- 
correlated RMF model. The diagrammatic arguments of Gornyi et al. |^6| do not fulfill 
the purpose of such a check: although their calculation of the current-current correlation 
function seems valid, their calculation of the conductivity relies on a relation between the 
conductivity and the current- current correlation function that seems questionable for the 
delta-correlated RMF model, since it depends on the quantity r that formally vanishes 
within the SCBA. Instead their arguments would seem to be applicable when the RMF is a 
weak perturbation to sufficiently strong, short-range disorder, in which case a well-defined 
t does exist. A direct diagrammatic calculation of the first logarithmic correction to the 
conductivity, which is not simple even for short-range disorder, therefore remains to be 
attempted for the delta-correlated RMF model. 

We remark also that further analysis of the free energy, (|75|) , with its new term, is 
deserved. For example, the influence of the new term on the scaling of the conductivity at 
small conductance, and its interplay with the usual topological term, remains an interesting 
open question. The general expression ( |69| ) for the free energy of a disordered conductor in 
the presence of equilibrium currents also presents new possibilities for further investigation, 
possibly in different dimensionalities. 
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The method of derivation of the ballistic a-model in this paper, by which a controlled 
expansion of the free energy becomes possible through the smallness of the parameter r/r tr , 
opens a promising avenue for further development of ballistic field-theories. It is possible 
that progress in this direction may help to resolve the problem of "repetitions" |58| in the 
ballistic a-model for chaotic systems. 
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the course of this work and gratefully acknowledge the financial support of the 'Sonder- 
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APPENDIX A: NEGLECT OF A 2 TERM 

In this appendix we justify more carefully the neglect of the A 2 term in the original 
Lagrangian, Eq. ([8]). In fact this approximation is standard and has been used in previous 
analytic works on the subject |4T|j43|| . Indeed our approach applies also to the case of long- 
range impurities (scalar disorder) for which AMike terms are completely absent. 

Even so we investigate here more carefully the influence of the A 2 term. In fact its neglect 
was made only for convenience and is not necessary in the early stages of the calculation: 
instead, it may be included exactly in the disorder average. In doing so, C int of Eq. (|j) is 
replaced by 
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hit 



where 



Jr,r' 



where IL is a unit matrix in zj-space. The decoupling of the interaction term, £i nt , by a 
Q-matrix follows as before, as does the integration over the V'-nelds. The Lagrangian of 
Eq. fllOD is now replaced by 
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We see that an additional term has arisen inside the logarithm (second line of Eq. ( |A1| )). 
This term is second-order in the scattering amplitude, V y ', and acts to renormalize the 
chemical potential. Indeed, due to the presence of the Str operating on this term, we see 
this term has no effect on the saddle-point Q p = A p , since the latter is supersymmetric. 

Consider now fluctuations around the saddle-point, of the form parametrized by Eq. flIE|) 
(such that Qp = 1). The unusual feature of the additional term inside the logarithm of 
Eq. (|A1]) is that, due to the presence of the Str operator, it commutes with the rotation 
matrices U and V n . Performing an expansion of the logarithm in this term will therefore 
not give any additional contribution to the free energy to any perturbative order. Instead, 
the effect of the additional term is comparable to that of the (a)-modes, corresponding to 
fluctuations in Q p , and hence may be neglected consistently. 

An alternative argument for the neglect of the A 2 term may be formulated as follows: 
we may treat the A 2 term in the Lagrangian (|§p perturbatively and check that it gives only 
a parametrically small correction to the previously-derived terms in the free energy. 

A leading-order contribution of the A 2 term, that is fourth-order in A, may be estimated 
via a third-order cumulant expansion: this gives rise to the contribution to the Lagrangian 
of 

[ drdr'dr" (U(r) — A(r).V r r 3 ^(r)) U(r')— A(r') V r ,r 3 ^(r')) x 
J W mc / V mc / 

x U{v")— 2 A 2 {v'')^{v'')\) 
\ mc z J I 

Choosing a typical pairing of the A and if) fields gives an estimate of a typical resulting 
contribution to the free energy: 

6F 1 ~±[dr ^(n a -n 1 )^( n ,- ni) niniStr(g B1 (r)Q ltt (r)g ila (r)) 

r 111,112,113 



39 



In fact due to the combination of three Q's in the structure of 5F\, this terms does not 



contribute at all to the free energy of the form of Eq. (27), upon expanding in P m . In 
addition we note that it is second-order in the scattering amplitude V 1 ^: for this reason it 
may be checked that in any case it appears with a small factor (epTtj-)^ 1 <C 1 with respect 
to the usual collision term (see e.g. Eq. (|2T|)). 

The other leading-order contribution, that is fourth-order in A, appears from the second- 
order cumulant: 



5Co ~ / drdr' 



^(r)^ji4 a (r)^(r)) U{v')^- 2 A 2 {v')i,{v')\ V 
mc z J \ mc z J I 



Choosing a typical pairing as before gives an estimate of a typical resulting contribution to 
the free energy as 

5F 2 ~ 4" / dr *m-n 2 Str ( ( 5n 1 (r)Q n2 (r)), 

r nj ,n 2 

^-n 2 =P% J d 2 p3V£V£ +pF{xll _ n2) . 

We notice that the term 8F2 is again second-order in the scattering amplitude V % \ as com- 
pared to the usual collision term in the ballistic a- model (Eq. (pT|)) which is only first-order 
in V y *. As a result it is straightforward to check that the term 8F2 produces a renormal- 
ization of the collision term by a parametrically small factor (e^r tl .) _1 1. Higher-order 
cumulants of the A 2 term in the original Lagrangian will clearly produce only smaller renor- 
malizations of the collision term corresponding to higher powers of this small parameter. In 
this way we conclude that it is justified to neglect the A 2 term as claimed in the main text. 
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